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The 2///3/3-decay running sums for ^"Ge and ^^°Nd nuclei are calculated within a QRPA approach 
with account for deformation. A realistic nucleon-nucleon residual interaction based on the Brueck- 
ner G matrix (for the Bonn CD force) is used. The influence of different model parameters on the 
functional behavior of the running sums is studied. It is found that the parameter Qpp renormal- 
izing the G matrix in the QRPA particle-particle channel is responsible for a qualitative change in 
behavior of the running sums at higher excitation energies. For realistic values of Qpp a significant 
negative contribution to the total 2z^/3/3-decay matrix element is found to come from the energy 
' region of the giant Gamow- Teller resonance. This behavior agrees with the results of other authors. 

(N ■ 

5_l ■ PACS numbers: 21.60.-n, 21.60.Jz, 23.40.-s, 23.40.Hc, 

Neutrinoless double beta (Oj^/3/3) decay, if observed, can allow one to extract the absolute neutrino mass from the 
f — ' decay rate provided that the corresponding nuclear matrix elements can be calculated reliably (see, e.g., 0,13). But up 
to now there has been no complete agreement among different many body methods calculating these matrix elements. 
Thus, further tests of the reliability of the calculated nuclear wave functions are needed. A probe that is of particular 
relevance is two-neutrino double beta (2t^/3/3) decay, as it connects the same initial and final nuclear ground states 
that 0vl3p decay does. The corresponding nuclear transition amplitudes M^!^p have been determined accurately in 
1— H ■ several nuclei of interest (see, e.g. [3| and references therein), 
y , The total 2i//3^-decay matrix element AP'^ is dominated by the double Gamow- Teller (GT) transition matrix 
^ ' element Mq^j.. The latter is determined by the amplitudes of the single GT transitions connecting the initial and final 
I— 1[ ground states with all l"*" states in the intermediate nucleus. The GT strength distributions in the first and second 
legs of the decay provide a more detailed, differential test of the calculations (while M^" is an integral quantity). It 
[ is customary to represent the partial contribution of all l"*" states with the excitation energies < lo to the total 
. M^'' in terms of the so-called 2i^/3/3 running sums M^''(a;). The running sum allows one to test whether so-called 
' single-state dominance (SSD) is realized in calculations. The SSD hypothesis states that the single transition via the 
(•~^ , lowest 1+ state in the intermediate nucleus can account for the entire M'^'^ . 

' For the three double-/3-decay systems A = 100, 116, and 128, in which the ground state of the intermediate nucleus 
' has J"^ = 1+, the single (3 decay in the second leg and the electron capture for the first leg provide additional tests of 
■ the quality of the theoretical calculations Recently, the GT strengths relevant for 2i^/3f3 decay have been studied 
^\ [ experimentally by means of different charge- exchange reactions [^-[1], thus providing important detailed information 
on the GT strength distribution. These measurements have shown that in some nuclei the SSD is indeed realized. This 
^ means that the total matrix element of 2://3/3 decay for those nuclei can be reconstructed fully from the measurement 
results only. For all the nuclei it has been obtained that the dominating contribution to JVP'^ comes essentially from 
the low-lying GT strength and in such a case charge-exchange reaction studies allow the provision of an upper limit 
^ ; on M^". 

■ - ' The 2i//?/3 running sums calculated for ''^Ca within the nuclear shell model (ol-fl^ consistently show a substantial 
overshoot of the total Mg"^ at excitation energies about 7 MeV, which then is compensated by a negative contribution 
coming from the region of the giant GT resonance (GTR). The total M^':^ — 0.054 MeV~^ calculated in Ref. [l^l is 
in excellent agreement with the corresponding experimental value, M^!^^ = 0.05 ± 0.01 MeV""'^ Q. The same feature 
was also observed in the QRPA calculations for ^^Ge and ^°°Mo 0. In contrast, recent calculations of 2z//3/3 running 
sums by Moreno et al. [1^, systematically performed within the QRPA with a schematic residual interaction for 
many nuclei, do not reveal any destructive contribution coming from the GTR region. Puzzled by this difference, 
we study in this Brief Report the influence of different QRPA model parameters on the functional behavior of the 



- ' — \ 



running sums by utilizing a QRPA approach developed in Ref. [lj|. The approach allows one to use, for deformed 
nuclei, a realistic effective interaction based on the Brueckner G matrix derived from the nucleon-nucleon Bonn CD 
force. We present an analysis of the 2i//3/3 running sums calculated for ''^Ge and ^^°Nd within the QRPA approach 
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of Ref. [3] . We find that the parameter gpp renormalizing the strength of the G matrix interaction in the particle- 
particle channel is responsible for the change in sign of the contribution to the total M^rp coming from high-lying 
states. Thus, different choice of the gpp parameter in different QRPA calculations may lead to a different, constructive 
or destructive, contribution to the total 2z//3/3-decay matrix element coming from high-lying 1"*" states. Therefore, a 
realistic value of gpp should be chosen thoroughly within a QRPA approach before a conclusion regarding the validity 
of the SSD is drawn. For example, quenching of the axial-vector constant gA can play an important role in fixing 
gpp. For realistic values of gpp in the present calculation (neglecting the quenching effect), a significant negative 
contribution to the total 2;y/3/3-decay matrix element is found to come from the energy region of the GTR. The values 
of gpp adopted in Ref. [l3| are apparently smaller than in the present calculations, and this may be a reason for the 
aforementioned differences in the functional behavior of the calculated running sums. 

A detailed description of the QRPA approach employed here as well as the choice of the model parameters is given 
in Ref. Here, we give only the formulas that are important for the present calculation of 2i^/3/3 running sums 

The matrix element M'qj, within the QRPA is given in the intrinsic system by the following expression : 

,t2u _ {Q^\~PK\K+.rnf){K+,mf\K+,m,){K+,m,\rK\Qt) 

Here, K is the projection of the total angular momentum onto the nuclear symmetry axis, characterizing nuclear 
excitations in the intrinsic system, {Q'^\j3^\K^ ^mf) and (X'^, TOi|/3^|0^) are the transition matrix elements of the 

GT operator (an explicit representation of these matrix elements as well as of the overlap {K^ ,mf\K'^ ,mi) is 
given in Ref. As in Ref. [IJ], we shift the whole calculated QRPA energy spectrum in such a way as to have the 

first calculated l"*" state exactly at the corresponding experimental energy. In this case the energy denominator in 
Eq. ([l} acquires the form w/f^mim/ = (w^:,™/ +'^i<',rni —w^f^i J/2-|-a)j^+ , with 0)^^+ being the experimental energy 
of the first 1+ state relative to the mean ground-state energy of the initial and final nuclei, and {Eq. + Eo^)/2 and 
^K,mf{^K,mi) being the calculated QRPA energies relative to the initial (final) nucleus. All the calculated running 
sums in this work are represented in terms of Hi. 

The 2vj3j3 running sum Mq^{uj) is defined as the sum, Eq. ([T]), truncated by the condition uJK,mimf £ 
In the present work we adopt the model parameters of Ref. A realistic residual interaction based on the 

Brueckner G matrix derived from the nucleon-nucleon Bonn CD force is used in the calculation, along with a schematic 
separable one. The residual interaction is renormalized in both the particle-hole and the particle-particle channels by 
means of multiplicative factors gph and gpp , respectively (corresponding strength constants in the case of the schematic 
force are x ^^^d k). The particle- hole renormalization factor gph is determined by fitting the experimental position of 
the GTR in the intermediate nucleus. For both ^^Ge and ^^°Nd the corresponding values = 1.15 for the realistic 
force and x" — 3.73/^4° '' MeV for the phenomenological separable force were found in Ref. [1J|. The particle-particle 
renormalization factor gpp is chosen so as to fit the experimental value Mq^^ for a given nucleus. The corresponding 
fitted values of gpp and also found in Ref. along with the choices of the deformation parameter /32, are listed 
in Table n 



TABLE I: Values of the deformation parameter /32 for initial (final) nuclei adopted in the calculations along with the fitted 
values of the p-p strength parameters g^^ (for the realistic Bonn-CD force) and (for a phenomenological separable force). 



Nucleus 


P2 


9pp 


k" (MeV) 


'"Ge C^Se) 


0.0 (0.0) 
0.10 (0.16) 


0.94 
0.99 


0.087 
0.091 


i^"Nd (i^"Sm) 
Def. I 
Def. II 


0.0 (0.0) 
0.37 (0.23) 
0.24 (0.21) 


1.11 
0.78 
1.35 


0.051 
0.033 
0.053 



To study the dependence of the shape of the calculated 2i^/3/3 running sums on different parameters of the residual 
interaction, we have calculated Mq^^(w) not only for the values of (x°) and gpp (k°) listed in TablelH but also for 
the p-h and p-p strengths quenched by a factor of 2. The running sums M^rp[uj) calculated with different combinations 
of the quenched and unquenched strengths are plotted in Figs. [1] and [51 Each panel in the figures shows four types of 
calculated running sums corresponding to different combinations of quenched and unquenched strength parameters 
(columns are labeled by values of gph , and rows by values of gpp in the calculation) . 

As shown in Figs. Il|2| this is the p-p interaction parameter gpp (k) that is responsible for the qualitative change in 
the behavior of the calculated 21^/3(3 running sums. For the quenched strengths gpp = g^p/2 {k = k;°/2), all calculations 
reveal a monotonically growing M'^j.iui). The total calculated M'^j, is, however, rather strongly overestimated for 
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these quenched gpp. For the reahstic, unquenched, gpp = g^^ {n = k^), a substantial overshoot of the total Mg";^ at 
low excitation energies is compensated by a negative contribution coming from the region of the GTR at w wlO MeV. 
The absence of such a pronounced overshoot in the case of ^^"Nd^-^^^Sm 2i//?/3 beta decay calculated with set I of 
the deformation parameters /32 can again be attributed to a substantially smaller fitted value of g^^ (k'^) in this case, 
see Table m In this case the difference in (32 between the initial and the final nucleus is large, and this leads to a large 
overall suppression of the calculated Moip by a small BCS overlap factor 



'Ge(P2=0)-> Se(p2=0) '''Gt(^=0A)->Se(^=0.l6) 




CO (MeV) CO (MeV) CO (MeV) CO (MeV) 

FIG. 1: (Color online) Running sums MQy(a;) for ^^Ge— >-^®Se 2ul3/3 decay, for different choices of deformation and strength 
parameters. The solid black line represents the QRPA results obtained with a realistic residual interaction {G matrix from the 
Bonn CD force), while the dashed (blue) line corresponds to the calculation with a schematic separable force. The left and 
right panels correspond to different deformations of the initial and final nuclei in the calculation. Each panel consists of four 
sub-panels forming 2x2 matrix, whose columns are marked by the strength of gpp (or k) taken in the calculation [the full or 
a half of the fitted value Qpp (k*^) given in Table |l] and whose rows are labeled by the strength of gph (or x) in the calculation. 

We have calculated the 2j//3/3-decay running sums for ^^Ge and ^^°Nd nuclei within the QRPA approach of Ref. [3| 
which takes nuclear deformation into account and implements a realistic nucleon-nucleon residual interaction based on 
the Brueckner G matrix (for the Bonn CD force). We have studied the influence of different QRPA model parameters 
on the functional behavior of the running sums within the QRPA approach. It has been found that the parameter 
gpp renormalizing the G matrix in the particle-particle channel is responsible for a qualitative change in behavior of 
the running sums. Therefore, a different choice of the gpp parameter in different QRPA calculations may lead to a 
different, constructive or destructive, contribution to the total 2i//3/3-decay matrix element coming from high-lying 1^ 
states. For realistic values of gpp a significant negative contribution to the total 2z^/3/3-decay matrix element is found 
in the present calculation to come from the energy region of the giant GT resonance. This observation is in accord 
with other results [I, l9l-[l2| . The values of gpp adopted in Ref. are apparently smaller than those in the present 
calculations, and this may be a reason for the different functional behavior of the running sums calculated in Ref. [l3| . 
Thus, a realistic value of gpp should be chosen thoroughly within a QRPA approach to draw a conclusion regarding 
the validity of the SSD. 
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